Thermoelectric three-terminal hopping transport through one-dimensional 
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A two-site nanostructure (e.g, a "molecule" ) bridging two conducting leads and connected to a 
phonon bath is considered. The two relevant levels closest to the Fermi energy are connected each 
to its lead. The leads have slightly different temperatures and chemical potentials and the nanos- 
tructure is also coupled to a thermal (third) phonon bath. The 3x3 linear transport ("Onsager") 
matrix is evaluated, along with the ensuing new figure of merit, and found to be very favorable for 
thermoelectric energy conversion. 

PACS numbers: 84.60.Rb, 72.20.Pa, 72.20.Ee 



I. INTRODUCTION 

In thermoelectric transport temperature differences 
can be converted to (or generated by) electric voltages. 
Such phenomena have already found several useful ap- 
plications. Current research is motivated by the need 
for higher performance thermoelectrics as well as the 
pursuit of understanding of various relevant microscopic 
processes (especially the inelastic ones). Theory^— pre- 
dicts that high values of the thermopower follow when 
the carriers' conductivity depends strongly on energy. 
Indeed, in bulk systems, the thermoelectric effects ne- 
cessitate electron-hole asymmetry, which is often rather 
small. However, in nanosystems, such asymmetry can 
arise in individual samples in ensembles with electron- 
hole symmetry on average. Moreover, inelastic processes 
and interference effects may play nontrivial roles in ther- 
moelectric transport £ It is known that the thermoelec- 
tric performance is governed by the dimensionless figure 
of merit ZT where T is the common temperature of 
the system and Z = aS 2 /(n e + K p h), with a being the 
electrical conductivity, S the Seebeck coefficient, and n e 
and Kph the electronic and the phononic heat conduc- 
tivities, respectively. Both n e and K p u can be smaller 
in nanosystems^ than in bulk ones, opening a route for 
better thermoelectrics. 

Mahan and Sofo 6 - have argued that the best thermo- 
electric efficiency can be achieved in systems where i) the 
energy width of the main conducting channel is very nar- 
row, and ii) the phonon thermal conductivity is as small 
as possible. It was suggested that ii) can be also real- 
ized in nanoscale composite structures where phonons are 
scattered by large variations in geometry and abundant 
interfaces^ while i) leads 6 to a very small K e /(S 2 <r), 
as indeed has been confirmed in studies of quantum dot 
arrays. 9 

Here we consider three-terminal thermoelectric trans- 
port in small one-dimensional (ID) nanosystems accom- 
plished via inelastic phonon-assisted hopping, and show 
that such processes lead to several nontrivial proper- 



ties. Although thermal transport properties of meso- 
scopic structures have been studied in the past^^ ' 10 ! 11 
investigations of three-terminal thermoelectric transport 
are just at their infancy^ Our main conclusions are 
drawn from the simple, but important, two localized- 
state junction in which hopping is nearest-neighbor. 13 
Later, we briefly discuss larger ID systems, which exhibit 
rather surprising features of the thermoelectric transport. 
We show that such systems may have a high figure of 
merit, as K e can become extremely small, while the ther- 
mopower, S, remains finite. Some of the thermoelectric 
transport coefficients we find correspond to transferring 
electric/thermal current via temperature difference be- 
tween the electron system and a suitable phonon bath. 
Hopefully, such systems can be achieved within current 
technology and be useful in applications. 

There are several related ideas in the literature. Ref- 
erence [3 presented an early, ingenious, way to cool a fi- 
nite 2D electron gas (which plays the role of the thermal 
bath) at low temperatures by elastic electron transitions 
to/from the leads. All the energies involved are only of 
order k^,T. Reference [ll] provided an experimental real- 
ization of some of the suggestions of Ref. [lj, with further 
analysis. Reference [IH demonstrated a quantum ratchet, 
converting the nonequilibrium noise of a nearby quan- 
tum point contact to dc current. Reference [TtI suggested 
a sophisticated Carbon nanotube structure, designed to 
extract energy from a discrete local oscillator at ultralow 
temperatures. The present work considers the full three- 
terminal case, where the energies involved can be larger 
than k^T, and a real reservoir can be cooled, not just 
one or several degrees of freedom. 



II. MODEL SYSTEM 

A. Hamiltonian 

The Hamiltonian, H = H e + H e - P h + H v h, consists 
of the electronic and phononic parts and the electron- 
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phonon interaction. The electronic part is (electronic 
operators are denoted by c and c') 

H e = ^2 E i c l c i + X/ e k(p) c l(p) c Hp) 

i k(p) 

+ ( S J 4 ,fc(p) c ! c fc( P ) + Yl J i,i+i4 c i+i + H - c -) • co 

i,k(p) i 

Here i labels the localized states, of energies -Ei (includ- 
ing usual Coulomb-blockade effects; i.e., it is assumed 
implicitly that a large Hubbard interaction confines the 
occupation of each level to be or 1) and k (p) marks 
the extended states in the left (right) lead, of energies 
(e p ) (all energies are measured from the common chemi- 
cal potential) . The matrix element coupling the localized 
states to each other is Ji_j, and those coupling them to 
the lead states are Ji,k( P )- All are exponentially decaying, 
with a localization length £, e.g., 



and Nb is the Bose function 



J.. 



i,k{p) 



a e exp 



J L(R) 



(2) 



with Xi and x^^ being the coordinates of the center of 
the localized states and the left (right) boundary, and the 
prefactor a e yielding the coupling energy. The electron- 
phonon interaction is 

H e-ph = X M V3 C l C j (°q + fl -q) + H ' C - ' ( 3 ) 

q 

where the phonon modes, of wave vector q and fre- 
quency ujq, are described by the operators a q , a q . Their 
Hamiltonian is H p h = w q a q a q (we use units where 
h = 1). The electron-phonon coupling is M^ij = 
c*e-ph exp(— \xi — Xj\/^), with a e - p h being the electron- 
phonon coupling energy. The transport through the sys- 
tem is governed by hopping when the temperature is 
above a crossover temperature, T x , estimated below in 
Sec. Ill CI for the most important two-site case. At lower 
temperatures the dominant transport is via tunneling. 
The two-site example of our system is depicted in the 
upper panel of Fig. [TJ 



B. Hopping and interface resistors 

The system described above bridges two electronic 
leads, held at slightly different temperatures and chemi- 
cal potentials, Tl, /il, and T R , /j R , such that the common 
temperature is T = (Tl + T r )/2. The golden-rule transi- 
tion rate Ty, between two localized states, located at Xi 
and Xj and having energies Ei < < Ej^- necessitates 
the inelastic electron-phonon scattering ([3]), and reads 

r ij =2nT in f i (l-f j )N B (E ji ) , (4) 

where Eji = Ej — Ei, the carriers' local Fermi function 



fi = 



exp 



V k B T t ) 



(5) 



exp 



f^q_) 



(6) 



determined by T p h, the temperature of the local phonon 
bath (see Fig. [TJ. We assume that this phonon bath is 
strongly coupled to a thermal reservoir and is thermally 
isolated as much as possible from the leads, such that 
its temperature is determined by that reservoir. On the 
other hand, phonons in the leads are in good thermal 
contact with the electrons there and share the same tem- 
perature. These assumptions are further elaborated upon 
in Sec. iHTCl In Eq. Q, T m = |M q ,y ;\ 2 v P h{\E ij I), where 
v p h is the phonon density of states. The linear hopping 



conductance at long distances (|xj 



high energies (\Ei\, \Ej\ » feT) of the bond ij it 
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- P h\ 2v ph{\ E ij\)' l lij 1 



Vh = ex P 



\Ej\ + \E,\ + \E l: 
2fc n T 



(7) 



As opposed to Eq. (f4j, the tunneling conduction from, 
say, site i to the left lead can be accomplished by elas- 
tic tunneling processes with a transition rate TiL = 
JiLfi[l-fL(Ei)], where -f iL = 2n\J. lk \ 2 i> L (E. l ) and f L and 
vl are the Fermi distribution and density of states of the 
left lead. The corresponding linear interface conductance 



is then GiL 



^VLiE^kvT)- 1 exp[-2|x iL |/e 



\Ei\/ (hsT)]. This conductance (and the interface con- 
ductance at the right lead) will be assumed to be much 
larger than the hopping conductance between the two 
localized states. 



C. The two-site case 

The thermopower in the hopping regime has been 
discussed by ZvyaginJ^ The simplest example is that 
of a two-site system — 1,2) depicted in Fig. [JJ 
which describes, e.g. a diatomic molecule^ or a series- 
connected double quantum dot.— In such a case trans- 
port is accomplished by nearest-neighbor hopping. As 
site 1 (2) is in a good contact with left (right) lead, we 
may assume that the local chemical potential and tem- 
perature there are Hl(r) an d T L ^ R y The transport is 
dominated by the hopping from 1 to 2 when the tem- 
perature is higher than T x . This temperature is esti- 
mated from the requirement that the elastic tunneling 
conductance across the system, Gtun, is comparable to 
the hopping one. The former is given by the trans- 
mission Ei=i,2 r iL{E)r iR (E)/[(E - E,) 2 + (T lL (E) + 
Tm(E)) 2 /4:], where the tunneling rates are T iL ( R )(E) = 
^\Ji,k(p)\ 2v L{R){E). Since site 1 (2) is coupled mostly 
to the left (right) lead, we use their perturbation-theory 
mixtures, governed by the small parameter 3\ijEi\. At 
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For an electron transferred from left to right, the bath 
gives an energy — E\ (E 2 ) to the left (right) lead, and 
thus the phonons transfer the energy E 2 \ to the electrons. 
A net energy of E = [E\ + E 2 ]/2 is transferred from left 
to right. Hence, the net electronic energy current, Iq, 
and the heat current^! exchanged between the electrons 
and the phonons, Tq are 



EI 



and Iji 



E 21 I N . 



(11) 
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b Phonon bath T ph 

FIG. 1. (Color online) a. A two localized-state (i and j, 
gray points) system coupled to two leads, of temperatures Tl 
and Tr, and chemical potentials {il and fiR (with the choice 
Mi > f-R an d Tl > Tr). The phonon bath temperature is 
Tph- The localized states are coupled (dotted lines) to the 
continuum of states in the leads, and are also coupled (the 
wavy line) to the phonon bath; b. The effective resistors rep- 
resenting the system: The straight (blue) arrows indicate the 
net electronic currents and the wavy (brown) one the phonon 
heat current, with Gil, G2R, and G12 being the conductances 
of the tunneling and the hopping resistors, respectively. 



low temperatures and for |I?i|,|.E 2 | ^ r l7 r 2 , where 
Ti = TiL + r^., we estimate 



G tun ~ e 2 I?rVl%(OK(0)£ 2 - 2 exp(-— ^ 



(8) 



where W is the system length between the leads. The 
hopping conductance is given by Eqs. j7]) (with i,j — 
1,2). Comparing those two, with exponential accuracy, 
the elastic tunneling mechanism can be important 2 ^ only 
when 

lE^ + \E 2 \ + \E 12 \ 



exp 



< exp 



W 



F12 



2k B T 

-y m3 > exp(2W70 , (9) 
giving k B T x - (|£?i| + \E 2 \ + \E 12 \)£/[4:(W - \x 12 \)}. 



The linear-response transport coefficients are obtained 
by expanding Eqs. (fTU|) and (jTTJ) to first order in ST = 
Tl - T R , 8n = [i L - [i R , and AT = T ph - T, 




Li K L 3 
L 2 L 3 K pe 



( Sp/e \ 

ST/T 
\ AT/T J 



(12) 



where I e — eljq is the charge current. All transport 
coefficients in Eq. (If 2[) are given in terms of the linear 
hopping conductance G [given by Eqs. 0], 



Lj — —E , L 2 — —E 21 

e e 



e e 2 



G— , 



G 



L 3 — -=EE 21 , K pe — -^E 21 



(13) 



Note that L 2 , L 3 , and K pe are related to E 2 i, and Eq 
vanishes linearly with the latter. 

The transport coefficients L 2 (£3) correspond to, 
e.g., generating electronic current {energy current) via 
the temperature difference AT.— When reversed, this 
process performs as a refrigerator: Electric current 
pumps heat current away from the phononic system 
and cools it down. In analogy with the usual two- 
terminal thermopower S, here we use the three-terminal 
thermopower— of this process, 



TG 



e krtT 



(14) 



Note that S p of our model can be very large as the energy 
taken from the phonons per transferred electron can be 
several times k B T . 



B. Two-terminal figure of merit, for AT = 



III. THREE-TERMINAL THERMOELECTRIC 
LINEAR TRANSPORT 



A significant feature of our setup is that the electronic 
heat conductance can vanish while the thermopower 
stays finite 



A. Transport equations 

The electronic particle current through the system is, 
as in Eq. Q allowing for the temperature and chemical 
potential differences, 

i N = r ia - r 21 . (10) 



T 2 

J1 e Ji e q w 



s = 



TG 



E 



e k B T 



(15) 



According to Ref. |fj, the largest two-terminal figure of 
merit is achieved in systems with the smallest K e /o~S 2 
(provided that S stays finite). Here this ratio vanishes, 
and then Z is limited by K p h- The latter can be minute in 
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nanosystems.^ Moreover, it can be reduced by manipulat- 
ing phonon disorder and/or phonon-interface scattering 
(avoiding concomitantly drastic changes in the electronic 
system). Our system is then expected to possess a high 
figure of merit. 



C. Three-terminal figure of merit 

The three-terminal geometry suggests novel possibili- 
ties for thermoelectric applications. For example, when 
AT < and S\i > 0, the setup serves as a refrigerator of 
the local phonon system, whose efficiency is given by the 
rate of the heat pumped from the phonon system to the 
electrical work invested, 



n = P Q e /{iM . 



(16) 



Consider first the special situation with E = 0, where 
L\ = = L3 = 0. For a given AT, Sfi is adjusted to 
optimize the efficiency, yielding 



V : 



?7o(2 + ZT- 2VZT+1)/{ZT) 



(17) 



where 770 = T/|AT| is the Carnot efficiency, and the new 
figure of merit is 



ZT = L 2 2 /(GK p , 



(18) 



Inserting here Eqs. (113)) yields ZT — >■ 00 upon neglect- 
ing the "parasitic" conductances, discussed below. When 
i? / 0, such an optimization can be achieved by setting 
ST = 0. 

In reality, ZT must be finite. To calculate a more re- 
alistic efficiency, we generalize Eq. (|12|) by adding the 
elastic transmission, the tunneling conductance G e i, to 
the hopping conductance G, and the elastic components, 
Li,eZ and K® e[ , to L\ and K®. (The elastic transmission 
does not contribute to L 2 , T3, and K pe , which are related 
to the heat transfer between the electronic and phononic 
systems^) We also include the phonon heat conductance 
K p , replacing Iq by Iq, the total heat current from the 
left to the right lead (ST is now also the temperature 
difference for the phononic systems in the left and right 
leads). Due to the absence of phonon-drag effects in lo- 
calized electronic systems, the temperature difference ST 
should not contribute to other currents beside Iq. Fi- 
nally, there are phononic heat flows from the two leads 
to the system being cooled. Hence the numerator of Eq. 
(fTo) is replaced by Iq — K pp AT/T, where K pp describes 
the phononic heat conductance in such processes. 

Following the same procedure as above, the efficiency 
is optimized by adjusting <5/i at ST = 0. The result is 
similar, except that the figure of merit is modified, 



ZT 



7-2 

L 2 



(G + G el )(K pp + K pe )-Lj 



G 



G el K pp y< 



GK„ 



(19) 



This has a straightforward physical interpretation: The 
wasted work is due to the elastic conductance and the 
unwanted heat diffusion, and ZT is limited by the ratio 
of the waste to the useful powers. In nanosystems K pp 
can be limited by the contact between the system and the 
leads. Hence the ratios can be made small and ZT can 
still be large. The three-terminal device can also serve as 
a heater and as a thermoelectric battery, where the same 
figure of merit describes the efficiency^ 
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FIG. 2. (Color online) a. The conductance, ln(G/Go), and 
thermopower, eS/ks; the abscissa gives the number of compu- 
tations, iV com . The parameters are W = 800, Lm = 20, £ = 
3.3, in units average nearest-neighbor distance, To = 3000, 
T — 20, and the energy-band width is 1090 (To and Lm are 
the Mot1s-*» temperature and length, defined in the text), b. 
The relative change in S as a function of the system length 
1T/(4Lm), obtained by averaging over 10 6 random configura- 
tions. Parameters (except W) are the same as in a. 



D. Longer ID systems 

For a chain of localized states, the picture is similar 
though slightly more complex. Consider first nearest- 
neighbor hopping, where the system is a chain of resis- 
tors. The same considerations as in the two-site case 
[in which the energy transferred is determined by site 
1 (2) and the left (right) lead] hold here for the left- 
most, t, (rightmost, r) localized state and the left (right) 
lead. Hence, the thermoelectric transport is described 
by Eqs. d]!-®, with E x (E 2 ) replaced by E e (E r ). 
In particular, the thermopower coefficients S and S p are 
completely determined at the left and right boundaries, 
despite the fact that transport coefficients are usually 
determined by both the boundaries and the "bulk" . 

In the variable-range hopping regime, the result is simi- 
lar: S and S p are determined by the resistors closed to the 
left and right boundaries, within a distance comparable 
to the Mott length. This observation is confirmed by nu- 
merical simulations. In Fig. [5Ji we plot the thermopower 
S and the conductance G for different random configura- 
tions in which the energies Ei and locations Xi of the sites 
are random: The Ei are chosen from a uniform distribu- 



tion in the range [-£ roax , E msx ], with E n 



545 (units 



are defined in the figure caption) being larger than the 
hopping energy, ~ (TqT) 1 / 2 , determined by the Mott2S 
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temperature To which is of the order of the level spacing 
on scale £). The locations chosen from a uniform 

distribution in the range [0, W], where W is the length of 
the sample in units of the nearest-neighbor distance. The 
conductance of the whole network is calculated by solv- 
ing the Kirchhoff's equations, 22 which yields the currents 
through each bond for a given source-drain bias. The 
thermopower S is obtained from the particle current, In, 
and the heat current Iq via the relation S = Iq/ '{I^eT). 

The heat current Iq (particle current) is calculated 
by summing over the heat currents (particle currents) 
flowing through all the bonds connected with the leads. 
In each such bond (iL) [(iR)] which connects the i— th 
localized state to the left (right) lead, the heat current 
flow is I™ = EiI W {I M = with ^ 

being the particle current in that bond. The total heat 
and particle currents are Iq = 0.5 ^(Tq^ + Iq R ^) and 

!n = J2t W = J2t 4? ' respsectively. 

At each N conl — th computation, with N com being an 
odd number, a new random resistor network is generated. 
At the subsequent, -/V com + 1— th , computatio only the 
middle half part [W/4, 3W/4] of the network is replaced 
by a new random configuration, while the parts close to 
the left and right boundaries, [0, W/4] and [3W/4.W], 
are not modified. It is seen from Fig. [2^ that the con- 
ductance G changes dramatically when the central part 
is modified, whereas the thermopower is practically im- 
mune to modifications of the central part. To further 
study the sensitivity of the thermopower to the sites that 
are a distance larger than W/4 away from the two inter- 
faces as a function of W, we plot the relative change of the 
thermopower S as a function of W in Fig. [2b- This rela- 
tive change is defined as \S 2n +i-S2n+2\/\S2n+i + S2n+2\, 
where S^n+i an{ i S 2n +2 denote the thermopowers cal- 
culated in the In + 1 and 2n + 2-th computation, re- 
spectively. The results are obtained by averaging over 
10 6 random configurations. It is seen that the relative 
change in S decays exponentially with increasing W, im- 
plying that sites located several Mott hopping distances 
Lm away from the boundaries have a negligible effect on 
the thermopower S. The Mott length, Lm, is of the or- 



der of [^/(vkBT)} 1 ^ d+1 \ where v is the density of states 
at the Fermi level. 



IV. SUMMARY AND CONCLUSIONS 

In this work we studied the three-terminal thermo- 
electric transport and thermopower mainly in simple 
ID hopping systems in the linear-reponse regime. We 
worked out the figure of merit for the three-terminal 
thermopower and expressed it as a function of the three- 
terminal thermoelectric transport coefficients. We ob- 
tained expressions for the thermoelectric transport coef- 
ficients in the simple two-site case. The system studied 
exhibits a large thermopower and high figure of merit in 
the appropriate cases. We analyze the conditions for high 
figure of merit in reality. 

For longer ID chains, we found that, contrary to in- 
tuition based on the usual conductances, the thermo- 
electric transport coefficients in hopping systems are 
solely determined by the states which are close to the 
interfaces (approximately, within the relevant hopping 
length) . More details on these surprising results and their 
generalizations will be given elsewhere. Finally, it should 
be emphasized that all the results obtained here are in 
agreement with the systematic microscopic derivations 
which will appear in consequent work. Both Eq. (|8]) and 
Eq. (|13p agree with the results derived from the non- 
equilibrium Green function method. 
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